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The discovery of dynamic memory effects in the magnetization decays of spin glasses in 1983 
marked a turning point in the study of the highly disordered spin glass state. Detailed studies of 
the memory effects have led to much progress in understanding the qualitative features of the phase 
space. Even so, the exact nature of the magnetization decay functions have remained elusive, causing 
confusion. In this letter, we report strong evidence that the Thermoremanent Magnetization (TRM) 
decays scale with the waiting time, tw By employing a series of cooling protocols, we demonstrate 
that the rate at which the sample is cooled to the measuring temperature plays a major role in the 
determination of scaling. As the effective cooling time, t%^^ , decreases, j- scaling improves and for 

tl^f < 20s we find almost perfect j- scaling, i.e full aging. 
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PACS numbers: 75.50 Lk 

Since the discovery of aging effects in spin glasses ap- 
proximately twenty years agO|lj||^], much effort has gone 
into determining the exact time dependence of the mem- 
ory decay functions. In particular, memory effects show 
up in the Thermoremanent Magnetization (TRM) (or 
complementary Zero-Field Cooled (ZFC) magnetization) 
where the sample is cooled through its spin glass tran- 
sition temperature in a small magnetic field (zero field) 
and held in that particular field and temperature con- 
figuration for a waiting time, At time t^, a change 
in the magnetic field produces a very long time decay 
in the magnetization. The decay is dependent on the 
waiting time. Hence, the system has a memory of the 
time it spent in the magnetic field. A rather persua- 
sive argument suggests that for systems with infinite 
equilibration times, the decays must scale with the only 
relevant time scale in the experiment, i^,. This would im- 
ply that plotting the magnetization on a i/tw axis would 
collapse the different waiting time curves onto each other. 
This effect has not been observed. 

What has been observed Q is that the experimentally 
determined magnetization decays will scale with a mod- 
ified waiting time, (t^,)''. Where is a fitting parame- 
ter. For /X < 1 the system is said to have subaged. A 
/J, > 1 is called superaging and /J. = I corresponds to full 
aging. For TRM experiments a /i of approximately 0.9 
is found for different types of spin glasses 0|, over a 
wide range of reduced temperatures indicating subaging. 
At very low temperatures and temperatures approach- 
ing the transition temperature, /i is observed to decrease 
from the usual 0.9 value. Superaging has been observed 
in Monte Carlo simulations of spin glasses 0. This has 
led to confusion as to the exact nature of scaling. 



Zotev et al. I?! have suggested that the departures from 
full J- scaling, observed in aging experiments, are mainly 
due to cooling effects. In a real experimental environ- 
ment, the situation is complicated by the time it takes 
for the sample to cool to its measuring temperature. An 
effect due to the cooling rate at which the sample tem- 
perature approaches the measuring temperature has been 
known 8, 9]. This effect is not trivial, it does not con- 
tribute a constant time to tw 

Another possible explanation for the deviation from 
full aging comes from the widely held belief that the 
magnetization decay is an additive combination of a sta- 
tionary term (M pf„,t. = A{To/t)") and an aging term 
(M = /(j^) )lia,ElE3- Subtraction of a stationary 
term, where tq is a microscopic spin flipping time, A is a 
dimensionless constant and a is a parameter determined 
from x" measurements, was shown to increase fi from 0.9 
to 0.97[l2. 

In this letter we analyze effects of the cooling time 
through a series of different cooling protocols and we 
present the first clear and unambiguous experimental ev- 
idence that the TRM decays scale as -^{i-e. full aging). 

Three different methods have been regularly employed 
to understand the scaling of the TRM decays. The first 
and simplest is to scale the time axis of the magneti- 
zation decay with the time the sample has spent in a 
magnetic field (i.e.j^)^]. If the decays scale as a func- 
tion of waiting time it would be expected that the decay 
curves would overlap. This has not yet been observed. 

A second more sophisticated method was initially de- 
veloped by Struik 13] for scaling the dynamic mechani- 
cal response in glassy polymers and first applied to spin 
glasses by Ocio et al. 15] . This method plots the log of the 
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FIG. 1: The cooling time(t') start when Tg is crossed. For CPl-2 the temperature is never dropped below Tm- On the contrary 
for CP5-6 the temperature is not allowed to go above Tm- The effective cooling time, t^-^^ is associated with a peak in the 
ZTRM S(t). The S(t) for each each cooling protocol are plotted with the the same y-scale. 



reduced magnetization M/AIfc {Mfc is the field cooled 
magnetization), against an effective waiting time ^ = p- 
where 



A=^[(l + A)l-._l]; ^<1 



or 



X = t^ log[l + — ]; 



(1) 



(2) 



A value oi fi = I would correspond to perfect t/t^ scaling. 
Previous values of fx obtained on the decays using this 
method have varied from .7 to .940 for a temperature 
range .2 < < .95. A value of /i < 1 is called subaging. 
Finally a peak in the function S(t) 



S{t) = 



1 dM(t) 
H d[Log,o{t)] 



(3) 



as a function of time has been shown to be an approx- 
imately linear function of the waiting time 2] . This peak 
occurs at a time slightly larger then the waiting time, 
again suggesting possible subaging. 

In this study we use all three of the above scaling pro- 
cedures to analyze the data we have produced with dif- 
ferent cooling protocols. All measurements in this letter 
were performed on our homebuilt DC SQUID magne- 
tometer with a Cu.94Afn.06 sample. The sample is well 



documented|l4l| and has been used in many other stud- 
ies. The measurements described in this letter were per- 
formed in the following manner: The sample was cooled, 
in a magnetic field of 20 G, from 35 K through its transi- 
tion temperature of 31.5 K to a measuring temperature of 
26 K. This corresponds to reduced temperature of .83 Tg. 
The sample was held at this temperature for a waiting 
time tu,, after which time the magnetic field was rapidly 
decreased to OG. The resulting magnetization decay is 
measured Is after field cutoff to a time greater than or 
equal to 5tw . The only parameters we have varied in this 
study are tw and the rate and profile at which we cool the 
sample through the transition temperature to the mea- 
suring temperature. The sample is located on the end of 
the sapphire rod and sits in the upper coil of a second 
order gradiometer configuration. The temperature mea- 
suring thermometer is located 12.5 cm above the sample. 
Heat is applied to the sample through a heater coil lo- 
cated on the same sapphire rod 17 cm above the sample. 
Sample cooling occurs by heat transfer with the He bath 
via a constant amount of He exchange gas, which was 
previously introduced into each chamber of the double 
vacuum jacket. We have measured the decay time of our 
field coil and find that we can quench the field in less then 
0.1ms. We have also determined that without a sample, 
our system has a small reproducible exponential decay, 
that decays to a constant value less than the system noise 
within 400 seconds. In order to accurately describe our 
data we subtract this system decay from all of the data. 

In this paper we present TRM data for eight waiting 
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FIG. 2: In the first column the TRM decays curves for ~ 
50s, 100s, 300s, 1000s, 3600s, 6310s, 10000s are plotted for the 
different cooling protocols. In the middle column the decays 
are scaled with the waiting. The y-axis for columns one and 
two are the same. The decays collapse onto each other as t^-^^ 
decreases. In the third column the decays are scaled with /j.. 
Where A from eq0one is used. 

times {U,= 50s, 100s, 300s, 630s, 1000s, 3600s, 6310s, and 
10000s). The same TRM experiments were performed for 
six different cooling protocols. In this paper we use four 
of the cooling protocols. Figure 1 (top row) is a plot 
of temperature vs. time for four of the cooling proto- 
cols. These different cooling protocols, were achieved, by 
varying applications of heat and by varying the amount 
of exchange gas in the vacuum jackets. A more detailed 
description of the cooling protocols will be given in a 
foUowup publication. In Figure 1 (bottom row) we plot 
S(t) (EqEl) of the ZTRM protocol (i.e. zero waiting time 
TRM) in order to characterize a time associated with the 
cooling protocol, t^-^-^ ■ As observed in Figure 1 we have 
achieved effective cooling times ranging from 406s down 
to 19 seconds. These times can be compared with com- 
mercial magnetometers which have cooling times in the 
range of 100-400s. 

In Figure 2, we plot the data for the TRM decays (first 
column) with the four cooling protocols. It should be 
noted that the magnetization (y-axis) is scaled by the 
field cooled magnetization. 

The second column is the same data, as column one, 
with the time axis (x- axis) normalized by tw It can be 
observed for scaling (column 2) that as the effective 
cooling time decreases the spread in the decays decreases 
giving almost perfect scaling for the 19 second cooling 
protocol. 

The last column in Figure 2, is the data scaled, using fi 



scaling which has previously been described. It has long 
been known that the rate of cooling affected fi scaling and 
that fi scaling is only valid in the limit t^, » t^-^^ . We 
find this to be true and that the limit is much more rig- 
orous than previously believed. To determine /i scaling, 
we focused on applying this scaling to the longest waiting 
time data (i.e., = 3600s, 6310s and 10,000s). For the 
largest effective cooling time data, t^^^ = 406s, we find 
that we can fit the longest waiting time data with a fi 
value of .88. This is consistent with previously reported 
values of We do find, however, that TRM data with 
waiting times less then 3600s do not fit on the scaling 
curve. We find that scaling of the three longest wait- 
ing time decays produces n values which increase as t^-^^ 
decreases. We also find that as f^^^ decreases, the data 
with shorter t^j begins to fit to the scaling better. It can 
be observed that at fjj^— 19s we obtain almost perfect 
scaling for all of the data with a value oi ji = .999. How- 
ever we find we can reasonably fit the data to a range of 
/i between .989-1.001. The fitting for the large decay 
curves, — 3600s, 6310s and 10,000s, is very very good. 
Small systematic deviations, as a function of f^, occur 
for < 3600s with the largest deviations for 1^^= 50s. 
Even with an effective cooling time two orders of magni- 
tude less then the waiting time, one sees deviations from 
perfect scaling. We have also scaled the data using Eq. 
2. We find no noticeable difference between the quality 
of this fit and the quality of the fi — .999, for t^^f = 19s, 
shown in Figure 2. Data with longer cooling times cannot 
be fit with Eq. 2. We therefore conclude that full aging 
is observed for the long t^; data using the t^^^ — 19s 
protocol. 

It is clear from Figure 1 (bottom row), that the effect 
of the cooling time has implications for the decay all the 
way up to the longest time measured, 10,000 seconds. 
The form of the S(t) of the ZTRM is very broad. The 
S(t) function is often thought of as corresponding to a 
distribution of time scales (or barrier heights) within a 
system which has infinite equilibration times (barriers). 
The peak in S(t) is generally associated with the time 
scale (barrier height) probed in time t^,. In Figure 1, 
we observe that for the larger effective cooling times, the 
waiting times correspond to points on or near the peak 
of S(t) for the ZTRM. We therefore believe that for the 
larger effective cooling times there is significant contam- 
ination from the cooling protocol over the entire region 
of t^s used in this paper. Only for t^^^ = 19s coohng 
protocol do we find that the majority of tu,s occur far 
away from the peak in the S(t). 

All the data in figure 3 used the cooling protocol with 
fjf = 19s. In Figure 3a we plot the S(t)(ZTRM) for 
f^f^ — 19s with arrows to indicate the waiting times for 
the TRM measurements. It can be observed that after 
approximately 1000 seconds the slope of the S(t) function 
decreases, possibly approaching a horizontal curve, which 
would correspond to a pure logarithmic decay in M(t). If, 
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different types of scaling we have performed, with and 
without the subtraction of the weaker power law term. 

We find that even for t^-f-l" = 19s the peak in S(t) for 
> 1000s occur at a time larger then tu,(fig 3b). We 
find that we can fit the effective time associated with the 
peak in S(t) to fjf = tl^^. 

In summary, we have performed TRM decays over a 
wide range of waiting times (50s - 10,000s) for six differ- 
ent cooling protocols. We find that as the time associ- 
ated with the cooling time decreases, scaling of the TRM 
curves improves in the j- scaling regime and in the ^ 
scaling regime. In /j, scaling we find that as the effective 
cooling time decreases /j, increases approaching a value of 
.999 for fjf = 19s. For the t^Jf = 19s TRM decays, we 
find that subtraction of a small power law term (A=.06, 
a=.02) slightly improves the scaling. It is however likely 
that the small systematic deviations of the t^^-^ — 19s 
data as a function of tw are associated with the small 
but finite cooling rate. 

The authors would like to thank V. S. Zotev, E. Vin- 
cent and J. M. Hammann for very helpful discussions. 



FIG. 3: Tlie S(t)(ZTRM) for tc = 19s is shown with 
the different t^s marked(a). In figure b. the S(t) for all 
decays(t^-''^ = 19s) are scaled with tw The decay curves 
are scaled with tm(c) and /i-Scaled(e), compared to the same 
scaling for decays with a stationary part subtracted(d and f ) . 
For stationary part we use A = 0.06, a — 0.02 and tq = 10~^^ 

on the other hand, the slope is continuously changing this 
part of the decay may be described by a weak power law. 
Either way, this region would correspond to aging within 
a pure non-equilibrated state. We believe that the long 
waiting time data occurs outside the time regime that 
has been corrupted by the cooling time and that this is 
the reason that we have, for the first time, observed full 
aging. 

It has been suggested that subtraction of a station- 
ary coinponent of the magnetization decay will improve 
scaling 1 12j. The very long time magnetization decay is 
believed to consist of a stationary term that is thought 
to decay as a power law. We fit a power law, M{t) — 
A{ToltY to the long time decay (1000s-5000s), of the 
ZTRM for t^=19s. Using To = IQ^^'^s, we find a = .07 
and A — .27. Subtracting this power law form from 
the magnetization decay destroys scaling. We find that 
the subtraction of a much smaller power law term with 
A=.06 and a=.02 slightly improves scaling at both short 
and long times. While the a values for the two different 
power law terms we have fit to are quite different, both 
values fall within the range determined from the decay 
of x"E3- In Figure 3(c-d) and 3(e-f), we plot the two 
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